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Abstract 

We have developed a formulation of density functional perturbation theory for the calculation 
of vibrational frequencies in molecules and solids, which uses numerical atomic orbitals as a basis 
set for the electronic states. The (harmonic) dynamical matrix is extracted directly from the first 
order change in the density matrix with respect to infinitesimal atomic displacements from the 
equilibrium configuration. We have applied this method to study the vibrational properties of a 
number of hydrogen-related complexes and light impurities in silicon. The diagonalization of the 
dynamical matrix provides the vibrational modes and frequencies, including the local vibrational 
modes (LVMs) associated with the defects. In addition to tests on simple molecules, results for 
interstitial hydrogen, hydrogen dimers, vacancy-hydrogen and self-interstitial-hydrogen complexes, 
the boron-hydrogen pair, substitutional C, and several 0-related defects in c-Si are presented. The 
average error relative to experiment for the ~ 60 predicted LVMs is about 2% with most highly 
harmonic modes being extremely close and the more anharmonic ones within 5-6% of the measured 
values. 

PACS numbers: 63.20.Pw,71.15.-m,71.15.Mb,71.55.Cn 
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I. INTRODUCTION 



The knowledge of the structures of impurities and defects is an essential prerequisite for 
understanding the electrical and optical changes that these complexes induce in semicon- 
ductors such as crystalline silicon0'l. The presence of light impurities such as H, B, C or O, 
results in the appearance of infra-red (IR) or Raman active local vibrational modes (LVMs) 
usually well isolated from the frequency range of the phonons of the host material. The 
observation of LVMs coupled with isotope substitutions and uniaxial stress measurements 
provide precious information about the type and number of impurity atoms involved and the 
symmetry of the defect. However, these data are rarely sufficient to identify unambiguously 
the defect. Since the early days of Steini, a large number of vibrational modes have been 
identified through the interplay of experiment and theory. 

The calculation of LVMs at the ab-initio level provides a critical link between theory and 
experiment. This is particulary true in the case of hydrogen, since it binds covalently in the 
immediate vicinity of many impurities and defects thus giving rise to a number of LVMs in 
the range ~ 800 to ~2200 cm~^. Other common impurities in Si which produce LVMs are 
B, C, O, but any element lighter than Si can in principle be observed by LVM spectroscopy. 

The computation of systematically accurate vibrational frequencies is a challenge for 
first principles theory, given their sensitivity on details of bonding geometry and electronic 
structure. Typical accuracies in the calculated vibrational modes for light impurities in 
silicon in former works are within 3-10% of the experimental data, which means in some 
cases a deviation of over 100 cm~^. 

Various approaches have been used to calculate LVMs, from semiempirical modelsi, to 
ab-initio Hartree- Fockl and density functional theoryi'lH. In most cases, frequencies are 
calculated in the spirit of the frozen phonon approximation. One computes the total energy 
of the system in the equilibrium configuration (that in which the forces acting on the atoms 
are zero) and then for small displacements of selected atoms (either individually or in the 
direction of a normal mode, if this is known). The actual value of the atomic displacement 
(typically a few hundredths of an A) and the response of the nearby atoms are parameters 
chosen by the user. One can either fit the energy vs. displacement to a polynomial and 
extract a specific vibrational model'!, or compute the dynamical matrix by finite differences!. 
When a few specific modes is all that is needed, only the movement of the atoms involved 
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in those modes is considered. In these methods, it is not possible to completely isolate the 
harmonic contributions from the anharmonic ones, since finite displacements always involve 
some anharmonic effects. For this reason, the frequencies obtained in this approach are 
sometimes referred to as quasi-harmoni^. 

One can also calculate vibrational properties from constant-temperature molecular- 
dynamics simulations, for instance by extracting selected frequencies from the velocity- 
velocity autocorrelation functionS or by using more sofisticated spectral estimators, like the 
Multiple Signal Classification (MUSIC) algorithm.llllli This is computationally exhausting, 
since long molecular dynamics runs are required, but potentially very accurate^! This also 
allows the calculation of frequencies as a function of temperature. 

However, the calculation of vibrational frequencies does not neccesarily require the actual 
displacement of the atoms, as in the methods described above. Linear response theory (in 
particular through the application of perturbation theory in density functional theory) has 
been thoroughly used in the past0 to compute the response of the system to infinitesimal 
atomic displacements, and from that, the vibrational frequencies in the harmonic approxima- 
tion. This can be done with the only knowledge of the electronic solution in the equilibrium 
configuration. The advantage of this approach is that anharmonic effects are elliminated, 
and that no reference is needed to explicit finite atomic displacements. Besides, this ap- 
proach allows to compute phonons with arbitrary q vector in crystalline systems, without 
having to consider a supercell conmensurate with the periodicity of the phonon, as it is 
required in the frozen phonon and molecular dynamics approaches. 

We propose a method, based on density functional perturbation theory (DFPT), to com- 
pute vibrational frequencies in the harmonic approximation. We use a basis set of numerical 
atomic orbitals to expand the electronic wavef unctions, which makes the method computa- 
tionally very efficient, and allows us to calculate systems with a large number of atoms. We 
have applied it to make a systematic study of a number of defect centers in silicon, involving 
light impurities and their complexes with intrinsic defects (vacancies and self-interstitials). 
In most cases, the comparison of the calculated and measured vibrational frequencies is very 
favorable, improving on the results obtained by other approaches. 

The outline of this paper is as follows. We first discuss the theoretical method and the 
model used to describe the defects. Then we compare the vibrational properties obtained for 
a variety of complexes with experimental data as well as other first-principles calculations 
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in the literature. Finally, we discuss the results. 



II. METHODOLOGY 

A. Ground State Description 

In this work, we use the fully self-consistent ab-initio code SlESTAli'0. The electronic 
energy is obtained from density-functional theory (DFT)0lli within the local density approx- 
imation. The exchange-correlation potential is that of Ceperley-AlderEl as parameterized 
by Perdew and Zunger.0 Norm- conserving pseudopotentialS in the Kleinman-Bylander 
formal are used to remove the core electrons from the calculations. 

The valence electron wavefunctions are described with numerical linear combinations of 
atomic orbitals (LCAO) of the Sankey type,0 but generalized to be arbitrarily complete 
with the inclusion of multiple-zeta orbitals and polarization states.0 These orbitals are 
numerical solutions of the free atom with the appropriate pseudopotential, and are extrictly 
zero beyond some cutoff radius. This makes the calculation of the Kohn-Sham hamiltonian 
to scale linearly with the number of atoms,llllli allowing calculations in very large systems 
with a modest computational cost. 

In the present work, the basis sets include single-zeta (SZ), double-zeta (DZ), double- zeta 
plus polarization (DZP), and triple- zeta plus polarization (TZP). A DZP basis includes two 
sets of s and jo's plus one set of d's on Si, O, C or B, two s's and one set of p's on H. The 



radial cutoff of the atomic orbitals was determined as described in Ref. 0, with an energy 
shift of 0.5 Ry, and a split-norm of 0.15 for all the species except H, for which the split-norm 
was 0.5. The charge density is projected on a real space grid with an equivalent cutoff of 90 
to 150 Ry to calculate the exchange-correlation and Hartree potentials. The host crystal is 
represented by a periodic supercell of 64 host atoms, and the k-point sampling is reduced to 
the r point. This restriction appears to be quite sufficient for the calculation of vibrational 
spectra. Tests have been performed for selected defects in a 128 host atoms cell and the 
results are within a few wavenumbers of those obtained in the 64 host atoms cell. 

In order to determine the equilibrium structure of the defects studied, we have relaxed 
all the atomic coordinates with a conjugate gradient algorithm, reaching a tolerance in the 
forces of F^ax <0.01eV/A. The dynamical matrix for the whole cell is computed (see below) 
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from this ground state and its eigenfrequencies and eigenmodes obtained. 



B. Linear Response Theory 

A new implementation of DFPT has been developed to compute the electronic response 
to infinitesimal atomic displacements. As is well known from the "2n + 1" theorem in 
quantum mechanicsj^l the first-order change of the electronic wavefunction in a perturbative 
expansion allows the computation of the second-order change in the energies. This implies 
that only the properties of the unperturbed ground state are needed to obtain the linear 
response of the system. The extension of this theorem to DFT is that the knowledge of the 
first-order change in the electronic density determines variationally the second-order change 
in the energy. In this way, we obtain analytically the dynamical matrix from the gradient of 
the density relative to atomic displacement and, from it, the vibrational properties, without 
physically displacing any atom. 

We describe briefly here the key points of our formulation. A complete report will be 
published elsewhere. The change in the electronic wavefunction is obtained by solving the 
first-order perturbation expansion of the Schrodinger equation (Sternheimer equationB) 

SHi/jo^i + Ho5iJi = Seitpo^i + eo^iStpt (1) 

where -00,1 ^'I'e the ground state electronic wavefunctions and Sipi the first order perturbation 
of ■0j when an atom is displaced (if we consider atom a, this would be daipi)- As we expand 
these wave functions in terms of atomic orbitals: 

^o,(r) = Ec.M0M(r-R/.) (2) 
the derivatives can be directly written in terms of the derivatives of the atomic orbitals: 

dai^i{r) = Y^[daCi^(l)f,{r - R^) + Cif,^a4>^,{r - R^)] . (3) 

Only the orbitals (p^ centered on atom a appear in the last term which is equal to — V0^(r — 
R^). The change in the coefficients, daQ^, is obtained from equation |I|. di/ji is then used to 
compute the perturbation in the electronic density 

occ 

a„p(r) = Y.[daij;^i + ridai^i] . (4) 

i=l 



This allows the computation of the dynamical matrix, by explicit derivation of the forces on 



all the atoms (3 in the system (the expressions of which can be found in Refs. |T5| and |16] for 
our approach) with respect to the infinitesimal desplacement of one of them [a): 

(M„M,)V2D„, = = (5) 

We remark that only terms up to first order in the electronic wavefunctions appear in the 
resulting formulas. Note that only the linear effects are obtained in this method, which is 
consistent with the harmonic approximation implicitly assumed in the diagonalization of the 
dynamical matrix. Thus, one expects to obtain high-quality frequencies for the vibrational 
modes that are harmonic, but the frequencies of modes involving large anharmonic contribu- 
tions will be less accurate. Although phonons with arbitrary q vector can be obtained in the 
pressent approach, here we only calculate vibrations which are periodic with the simulation 
supercell {i.e., q = F), since we are interested in LVMs. 



III. RESULTS 



Tests of the method for free SiH4, CO, CO2 and H2 lead to very good agreement with 
experiment (see table 1). Using an appropriate basis set reveals to be essential to reproduce 
accurate frequencies. In most of the DZ set gives good values, but in some con- 

figurations a more complete basis is required. This is particularly true for bending modes. 
Oxygen likes to have polarization orbitals, and thus the frequencies are better when these 
are included. Note that the vibrational frequencies are more sensitive to the basis set size 
than the structural properties, such as bond lengths. In general, the largest improvements in 
the frequencies correspond to going from SZ to DZ, then DZ to DZP, but TZ basis produce 
only marginal improvements. 

A number of defects containing light impurities in c-Si are now considered. These are 
H at a bond-center (BC) site, H dimers (H2 and Hg), the hydrogenated vacancy (VIIn,n = 
1,2,3,4) and the saturated divacancy V2II6, the self-interstitial-hydrogen IH2 complex, the 
{B,H} pair, substitutional C, interstitial O, and two charge states of the A-center (oxygen- 
vacancy complex). These embrace a range of Si-X bonding configurations, with X=II, C, 
O and B. In most cases, we have tried different basis sets in order to check, or improve the 
accuracy of our calculations when the chemical properties are particularly complex. We show 
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TABLE I: Calculated and measuredc3 frequencies for free SiH4, CO, CO2 and H2 molecules with 
various basis sets. 





SZ 


SZP 


DZ DZP 


TZ 


TZP 


Expt. 


















2045 


2064 


2160 2153 


2173 


2159 


2191 


Ai 


1970 


1974 


2110 2116 


2131 


2125 


2187 


E 


800 


862 


921 929 


926 


929 


975 


T2 


701 


772 


806 818 


817 


821 


914 








H2 








A 


3619 


3670 


4194 4185 


4193 


4191 


4161 








CO 








A 


1681 


1823 


1885 2088 


1945 


2183 


2170 








CO2 








Ai 


2118 


2355 


2235 2277 


2224 


2394 


2349 


A2 


1107 


1216 


1200 1241 


1209 


1331 


1333 


E 


478 


558 


547 583 


560 


635 


667 



frequencies for the DZ and DZP basis sets. Although the diagonalization of the dynamical 
matrix gives all the (F point) modes in the cell, we present here only the stretching and some 
bending LVMs. We also obtain the eigenvectors which we use to determine the symmetry of 
the corresponding vibrational modes. Our results are compared with available experimental 
data, and with theoretical frequencies obtained by other authors using first-principle DFT. 

A. Hbc, H2 and 

At the BC site,i§ hydrogen exists in the +1 and the charge states. In the latter case, 
the odd electron does not participate in the bonding but resides in a non-bonding orbital 
primarily localized on the two Si atoms adjacent to the proton. Chemically, this is a 3-center, 
2-electron bond, very much like the type of H bonding occurring in boron hydrides. The 
bond is somewhat compressed because optimal relaxation (no Si second nearest neighbors) 
would likely result in a longer Si-H-Si bond and a frequency lower than observed. This 
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FIG. 1: Calculated structure of H2 complex in silicon. Dark spheres are Si atoms, white spheres 
are H. The two perpendicular arrows represent the Hab wag eigenmodes at 842 cm~^. The twofold 
degenerate wagging mode for H^c is found at 560 cm~^. 

suggests that as H moves along the trigonal axis, it tends to form something like Si...H-Si 
then Si-H-Si then Si-H...Si, a process which is highly anharmonic. The calculated frequency 
is in table ||. 

RamanS and IR^Ell measurements of H2 in silicon reveal a considerable softening of the 
stretching mode with respect to the frequency of H2 in the gas phase. A number of calcula- 
tions (for a review, see Ref. found the molecule to be stable at the tetrahedral interstitial 
(T) site. The electron affinity of the Si atoms surrounding H2 is at least partly responsible 
for a small charge transfer from H2 to its Si neighbors, which results in a weakening of 
the H-H bond. Even though the H-H stretch mode is not expected to be fully harmonic, 
our calculated frequency is close to the experimental one (table H). Note that the errors 
relative to experiment in the D2 and HD frequencies are very different than the error in 
the H2 frequency. This is also a clear feature of these frequencies when they are calculated 
dynamically from the v-v autocorrelation function.Ei 

The trigonal H2 defect! consists of one hydrogen atom close to the antibonding (AB) 
site, and the other near the BC site (see figure 1). The two H atoms are inequivalent. The 
Si-H^B bond length is slightly longer than the Si-H^c one (we obtain 1.580 and 1.510 A, 
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respectively) which gives rise to different stretch frequencies for the two atoms: 2135 cm~^ 
for B.BC and 1750 cm~^ for H^^. We also obtain two degenerate waggind modes, associated 
with the Hbc atom, at 560 cm~^, and two wagging modes with 839 and 843 cm~^, related 
to H^^. The latter should be degenerate, but small inaccuracies in the atomic relaxations 
render them at slightly different frequencies. 



TABLE II: Calculated and measured frequencies for Hbc (our calculation is spin polarized for 
HgQ and spin averaged for Hg^), H2 in the (100) alignment, and H2 = HbcHab- A DZ basis was 
used for all these complexes. The errors relative to experimental values are in parenthesis, (a) is 



Ref. gg, (6) is Ref. |3l|, (c) is Ref. g2| (d) is Ref. |3^, (e) is Ref. |3|, (/) is Ref. g5|, (g) is Ref. |36|, 



(h) is Ref. |3^, and (i) is Ref. 38. Were H2 a classical dumbbell, its wag modes would be at 731 
and 860 cm^^ (see discussion in Ref. 39). 
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B. Hydrogen and native defects 

A considerable number of IR and Raman lines are related to H-intrinsic defect complexes. 
It has been notecl@@ that vibrational modes above 2000 cm~^ are mainly related to H in 
vacancies, while those lines below 2000 cm~^ are predominant for the H-self-interstitial 
systems or H at AB sites. A large number of geometrical configurations may lead to very 
similar vibrational lines, making it difficult to identify these defects unambiguously. As noted 
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by other groupsaacSe^, the streching frequencies in VH^ (n = 1,2,3,4) systems increases 
with n due to the repulsive H-H interaction. Thus, the highest IR hne is that of VH4. The 
Si-H bonds point toward the center of the vacancy along the trigonal axes. 

In our calculations (table [TTTD , VH has monoclinic symmetry, and the H oscillates parallel 
to the (111) direction. In the orthorhombic VII2, the two equivalent H have stretching 
modes along the (100) and (001) directions. The frequencies for these modes are 2121 and 
2144 cm~^ respectively. VH3 has Cst, symmetry. The A singlet involves the movement of 
the three H atoms towards the vacancy, while in the twofold degenerate E mode one of the 
atoms moves in opposition. VH4 has symmetry. In addition to the threefold degenerate 
T2 mode at 2205, we obtain an IR-inactive singlet Ai mode at 2265 cm^^. 

The vibrational modes of V2H6 are almost identical to those of VH3: The fully saturated 
divacancy behaves very much like two weakly coupled VII3 complexes. The A2 singlet at 
2176 cm^^ induces a dipole along the (111) direction. In addition to this mode and the 
IR-active E doublet, we obtain two IR-inactive modes at 2186 cm~^ and 2134 cm~^. 

The IH2 complexi^ has two equivalent and weakly coupled hydrogen atoms, which yields 
two very similar stretching frequencies. Uniaxial stress measurements show that the two 
hydrogen atoms are equivalent. Our relaxed structure has almost symmetry, with the 
A mode higher than the B mode, confirming early results^. The deviation from perfect 
symmetry is due to the finite tolerance in the geometry optimization. This deviation is seen 
when comparing the IHD and IDH complexes: they should be identical but we find their 
frequencies to be off by 2 cm~^. Note that we reproduce the correct ordering for the bending 
modes of IH2 at 737 and 732 cm^^ (table 0). 



C. Oxygen,Carbon and Boron in Silicon 

Oxygen is a well-known impurity which is especially important in Czochralski-grown 
Si, and a considerable amount of effort was focused in understanding its properties, il We 
have computed the LVM frequencies for interstitial oxygen (Oj) and two charge states of 
the vacancy-oxygen complex (A-center). The results are in table 0. Frequencies for Oj 
were computed with the oxygen placed at the BC site, where the probability of finding 
this delocalized atom is maximun, and the classical harmonic potential can better describe 
the local modes. The IR-active A2U mode corresponds to the asymmetric-stretching mode. 
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while the Aig is the symmetric one. Eu mode involves the movement of nearest silicon atoms 
with no participation of the oxygenS. Finally, TablejV^ shows the triply degenerate mode 
of substitutional carbon as well as the LVMs associated with the {B, H} complex.l 

IV. CONCLUSIONS 

We have presented a new development of DFPT, using localized atomic wavefunctions as 
a basis set, and applied it to the study of LVMs for light impurities in silicon. In contrast to 
other methods, the dynamical matrix is computed analytically without actually displacing 
any atom from its equilibrium position. The calculations are based on the ground state 
density matrix as computed with the Siesta package. 

Tests of the method for free molecules (SiH4, H2, CO, and CO2) show that this approach 
is highly accurate in situations where the anharmonic contributions are small. Note that 
the frequencies are obtained at T = OK while experimental data are obtained at low, but 
non-zero, temperatures. 

We have used a variety of basis set sizes to describe the electronic wavefunctions. In 
most cases, a DZ basis is quite sufficient to obtain accurate atomic structures and vibrational 
frequencies. Larger basis sets such as DZP improve the frequencies in situations that involve 
more complex chemical bondings. The defects included here are Hbc? H2, Hg, VHn (with 
n=l,2,3,4), V2H6, IH2, the {B, H} pair, substitutional C, interstitial O and two charge states 
of the A center. These defects involve a wide range of bonding configurations. 

The average error of the 60 calculated modes relative to experiment is about 2%. In 
situations where large anharmonic contributions are present, the accuracy of the method 
decreases somewhat (5 — 6%). This occurs, for example, when H is close to a BC position. 
However, in most cases the calculated frequencies are in remarkable agreement (0 to 2%) 
with experimental data, implying that this perturbative approach is totally justified and 
that the ground state density matrix calculated with Siesta is very reliable. 
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TABLE III: Calculated and measured frequencies for stretching modes in VHn (n = 1, 2, 3, 4) and 
V2H6. (a) is Ref. (6) is Ref. and (c) is Ref. The error relative to experiment is in 
parenthesis. Our frequencies were obtained with a DZP basis. The measured values for the VH2D, 
VHD2 and VD3 as published in Ref. || are now believed to belong to a divacancy complexE§ and 
are therefore not listed here. 





expt." 


this work 


Ref.'' 




Expt."" 


this work 


Ref.^ 






VH 








VH4 




A' 


2038 


1971(-3%) 


2248(+10%) 


Ai 


2257= 


2265(+0%) 


2404(+6%) 






VD 




T2 


2222 


2205(-l%) 


2319(+4%) 


A' 


1507 


1418(-6%) 


1613(+7%) 














VH2 




Ai 


2250 


2251(+0%) 


2384(+6%) 


Ai 


2144 


2163(+1%) 


2316(+7%) 


E 


2224 


2205(-l%) 


2319(+4%) 


Bi 


2121 


2132(+1%) 


2267(+7%) 


Ai 


1620 


1594(-2%) 


1677(+3%) 






VHD 








VH2D2 




A' 


2134 


2135(-0%) 


2292(+7%) 


Ai 


2244 


2235(-0%) 


2364(+5%) 


A' 


1555 


1551(-0%) 


1641(+5%) 


Bi 


2225 


2204(-l%) 


2319(+4%) 






VD2 




Ai 


1628 


1603(-1%) 


1690(+4%) 


Ax 


1564 


1552(-1%) 


1658(+6%) 


B2 


1615 


1585(-2%) 


1663(+3%) 


Bi 


1547 


1532(-0%) 


1625(+5%) 






VHD3 








VH3 




Ai 


2236 


2221(-1%) 


2342(+5%) 


Ai 


2185 


2158(-1%) 


2318(+6%) 


Ai 


1636 


1613(-1%) 


1705(+4%) 


E 


2155 


2100(-2%) 


2256(+5%) 


E 


1616 


1584(-2%) 


1664(+4%) 






VH2D 








VD4 




A' 




2140 


2298 


Ai 


no-IR 


1623 


1721 


A" 




2101 


2256 


T2 


1617 


1584(-2%) 


1664(+3%) 


A' 




1520 


1632 














VHD2 








V2H6 




A' 




2121 


2277 


Ai 


2190^= 


2186(-0%) 




A' 




1534 


1646 


A2 


2191 


2176(-0%) 




A" 




1509 


1619 


E 


2166= 


2143(-1%) 








VD3 




E 


2165 


2134(-1%) 




Ai 




1547 


1661 T 
15 










E 




1510 


1619 











TABLE IV: Calculated frequencies with a DZ basis set compared with experimental and other 
theoretical results for IH2. (a) is Ref. ^ and (6) is Ref. 45. The error relative to experiment is in 
parenthesis. 





expt.'* 


this work other authors 


expt." 


this work 


other authors 








IH2 








ID2 


A 


1989 


2007(+l%) 2107''(+6%) 2145"( +8%) 




1448 


1440(-1%) 


1510^(+4%) 1540'^(+6%) 


B 


1986 


2004(+l' 


Vo) 2106'' (+6%) 2143" ( +8%) 




1446 


1438(-1%) 


1508^ (+4%) 1539" (+6%) 


B 


748 


737(-l' 


Vo) (A) 775'^ (+3%) 






609 


590" 


A 


743 


733(-l' 


Vo) (B) 768*^ (+3%) 






601 


583" 


A 




716 


{B) 736*^ 






566 


564" 


B 




711 


(A) 717'^ 






562 


555" 








IHD/IDH 










expt." 




this work 






other authors 




1988 




2005/2007(+l%) 






2106^(+6%) 


2144" (+8%) 




1447 




1440/1438(-1%) 






1509^ (+3%) 


1540" (+6%) 




746 




733/736(-2%) 








771" (+3%) 








714/714 








727" 
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TABLE V: Calculated and measured LVMs for interstitial O (Oj) and two charge states of the 
A-center (V0° and VO"). (a) is Ref. || for and Ref. || for VO^o/"), (b) is Ref. |5|, and (c) is 
Ref. For \0^^/^^ a DZP basis was used for O and its Si nearest-neighbors, and a DZ basis for 
the other Si atoms. For O,-, a DZP basis was used for all the atoms in the cell. 







expt.° 


this work 


Calc'' 


Calc^ 


















1136 


1131(-0%) 


1184(+4%) 1108(-2%) 


1098(-3%) 


0. 




618 


607(-2%) 


619(+0%) 621 (+0%) 


630(+2%) 




En 


518 


538(+4%) 


519(+0%) 518(+0%) 






Bi 


836 


861(+3%) 


839'' (+0%) 


843(+l%) 




Ai 


534 


546(+2%) 


548" (+3%) 


540(+l%) 




Bi 


885 


897(+l%) 


872" (-1%) 


850(-4%) 




Ai 


545 


558(+2%) 


532" (-2%) 


539(-l%) 



TABLE VL Calculated and measured frequencies for substitutional C and the {B, H} pair in Si. 



(a) is Ref. (6) is Ref. |8|, (c) is Ref. ||, (d) is Ref. |T|' (e) is Ref. |6|, and (/) is Ref. |6C. 
A DZP basis was used for C and its Si nearest-neighbors, and DZ for the other atoms. For the 
{B, H} pair, a DZ basis was used for all the atoms. 







expt. 


this work 


other authors 


c. 


Td 


607" 


631(+4%) 


662* (+9%) 




A 


1903^ 


1958(-F3%) 


1830'^(-4%) 


{B,H} 








1880'^(-1%) 




E 


652 


695(+6%) 
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